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ABSTRACT 

This  thesis  investigates  the  feasibility  of  simulating 
the  Landing  Signal  Officer  (LSO) /pilot  interaction  during 
the  approach  to  a  landing  on  an  aircraft  carrier.   A  simu- 
lator was  created  which  duplicated  the  LSO ' s  operational 
environment  through  the  use  of  computer-generated  visual 
displays.   The  LSO  and  the  pilot  were  placed  in  this  simu- 
lated carrier  approach  environment  by  1)  displaying  a 
representation  of  the  landing  area  plus  a  "meatball"  and 
angle  of  attack  information  to  the  pilot  while  2)  simul- 
taneously displaying  the  aircraft's  approach  to  the  LSO. 

Test  results  demonstrated  the  basic  feasibility  of 
simulating  the  LSO/pilot  interaction  and  its  application 
as  a  research  tool  in  studying  LSO  models,  wave-off  tech- 
niques and  landing  techniques . 
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I.   INTRODUCTION 

This  thesis  analyzes  the  feasibility  of  creating  a  simu- 
lation of  the  Landing  Signal  Officer  (LSO) /pilot  interaction 
during  an  approach  to  a  landing  on  an  aircraft  carrier. 

Many  aspects  of  the  carrier  landing  system  [Fig.  1]  have 
been  extensively  studied.   The  pilot-aircraft  interface  has 
been  investigated  with  respect  to  auto-pilot  and  Automatic 
Carrier  Landing  (ACL)  systems.   The  flight  dynamics  of  the 
aircraft  [Ref .  2]  and  the  dynamics  associated  with  the  carri- 
er have  also  been  well  documented.   Aircraft  disturbances 
resulting  from  wind  or  stack  gas  generated  burble  have  al- 
ready been  identified  and  reproduced.   Human  describing  func- 
tions have  replaced  the  pilot  in  some  simulators. 

The  one  area  lacking  quantitative  analysis  was  the  LSO/ 
pilot  interface  shown  as  a  path  parallel  to  the  pilot's  in- 
put cues  in  Figure  1.   The  first  systematic  study  of  the 
Landing  Signal  Officer  was  completed  by  Gail  J.  Borden  in 
1969  [Ref.  5,6].   This  study  stated  that  one-sixth  of  all 
carrier  accidents  reported  to  the  Naval  Safety  Center  during 
a  five-year  period  had  the  LSO  listed  as  either  a  contribut- 
ing or  causal  factor.   A  more  thorough  analysis  of  the  LSO ' s 
role  in  the  carrier  landing  system  was  suggested. 

Mr.  Borden  described  many  problem  areas  such  as  the 
difficulty  LSO's  encounter  in  judging  pilot/aircraft  approach 
performance  which  were  identified  by  means  of  data  received 


from  opinion  surveys.   He  suggested  that,  especially  in  the 
area  of  the  LSO ' s  preceptual  ability,  objective  validation 
of  the  opinion  data  was  necessary.   The  Naval  Air  Develop- 
ment Center  (NADC)  made  an  analytical  attempt  to  mathemati- 
cally simulate  the  behavioral  characteristics  of  the  LSO- 
pilot  combination  [Ref .  7] .   They  continued  with  this  program 
without  the  benefit  of  a  suitable  data  base. 

Mindful  of  its  potential  usefulness,  a  simulation  was 
created  and  subjectively  evaluated  to  see  if  the  LSO/pilot 
interaction  environment  could  be  reasonably  reproduced.   The 
physical  presence  of  the  pilot  and  the  LSO  would  avoid  the 
problem  of  trying  to  mathematically  model  them.   The  focus 
of  attention  was  primarily  on  displaying  to  the  LSO  a  picture 
of  the  approaching  aircraft  with  enough  realism  for  him  to 
detect  changes  in  lineup,  air  speed,  angle  of  attack,  and 
position  with  respect  to  the  glide  slope. 

Three  already  published  techniques  were  combined  in  the 
construction  of  the  interaction  simulation.   The  aircraft's 
equations  of  motion  derived  by  Etkin  [Ref.  2]  were  solved 
on  an  analog  computer.   Techniques  in  computer  processing  of 
three-dimensional  structures  developed  by  Desens  [Ref.  3] 
were  used  to  generate  the  displays  shown  to  the  LSO  and 
pilot.   The  general  procedures  used  by  Kahrs  [Ref.  4]  to 
display  a  carrier  landing  picture  to  the  pilot  were  used 
to  construct  the  pilot  side  of  the  simulation. 

The  interactions  of  the  final  combined  system  are  shown 
in  Figure  2.   The  Pilot,  with  the  visual  aid  of  a  computer 


generated  display  [Fig.  3]  consisting  of  a  runway,  angle  of 
attack  indicator,  and  a  "meatball",  maneuvered  the  aircraft 
during  the  approach  to  a  carrier  landing.   Concurrently, 
the  LSO  was  presented  with  a  visual  display  which  depicted 
the  aircraft's  approach  as  seen  from  the  LSO's  platform 

[Fig.  4].   An  analog  computer  solved  the  aircraft's  equa- 
tions of  motion  while  a  digital  computer  handled  the  spacial 
transformations  needed  for  the  graphical  display  [Fig.  2] . 
The  LSO's  input  to  the  pilot  was  aural. 

Prior  to  incorporation  into  the  simulation,  the  LSO's 
visual  presentation  was  tested  to  see  if  a  realistic  repre- 
sentation of  the  aircraft's  approach  could  be  displayed. 
An  "O.K."  pass  was  matched  for  realism  with  "O.K."  passes 
observed  and  photographed  from  the  aircraft  carrier  "Enter- 
prise".  A  test  program  was  then  set  up  without  the  pilot 
in  the  loop.   Each  pass  was  fixed  in  that  the  parameters 

(lineup,  angle  of  attack,  glide  slope,  or  air  speed)  were 
constant.   This  preliminary  program  was  then  tested  for 
validity  by  having  several  subjects  judge  various  approaches 
to  determine  parameter  sensitivity  and  real  world  discrep- 
ancies . 

The  computers  used  in  this  simulation  were  located  in 
Spanagel  Hall  on  the  campus  of  the  Naval  Postgraduate  School, 
Monterey,  California;  however,  the  techniques  used  were 
applicable  to  any  comparable  computer  systems.   The  computers 
used  were:   one  analog  (CI  5000);  one  32K  digital  (XDS  9300); 
and  two  graphics  systems  (AGT  10).  [Ref.  1]. 
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II.   ANALOG  AND  DIGITAL  COMPUTER  PROGRAMS 

A.   ANALOG 

Reference  2  contains  the  equations  of  motion  using  sta- 
bility axes,  small  disturbance  theory,  and  non-dimensionali- 
zation.   The  longitudinal  and  lateral  equations  of  motion 
were  reduced  to  a  form  useful  in  this  simulation.   These 
equations  and  the  assumptions  made  are  presented  in  Appendix 

A,  and  the  scaled  version  of  these  equations  as  solved  by 
the  CI-5000  analog  computer  is  shown  in  Appendix  B. 

Inputs  to  the  analog  computer  were  the  outputs  of  a 
stick  and  throttle  attached  to  a  seat  placed  in  front  of 
the  pilot's  graphics  displays.   The  outputs  of  the  analog 
program  were  the  inertial  coordinates  of  the  aircraft  (X', 
Y*,  Z'),  the  changes  in  Euler  angle  rotation  (8,  <j> ,  iJj)  and 
changes  in  the  Z'  component  of  velocity  (w) .   The  inertial 
position  of  the  aircraft  was  then  used  to  locate  the  air- 
craft * s   position  with  respect  to  the  LSO's  position.   The 
rotation  angles  were  used  to  properly  rotate  the  aircraft 
in  the  LSO's  display.   The  velocity  component  extracted  was 
used  to  define  the  airspeed  of  the  aircraft. 

B.  AXIS  MANIPULATION 

Using  the  seven  outputs  of  the  analog  computer,  the 
digital  computer  prepares  the  data  sent  to  the  graphics 
computer  for  display  [Fig.  6].   The  critical  factor  in 
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processing  the  inputs  for  graphical  display  was  the  rate 
at  which  the  picture  was  updated,  i.e.  loop  time.   Taking 
a  movie  projector  as  a  norm,  a  goal  of  not  more  than  18 
updates  per  second  was  set.   Therefore,  time-saving  de- 
vices were  incorporated  whenever  possible.   Each  display 
will  now  be  discussed  separately. 

A  thesis  by  R.  B.  Desens  [Ref .  3]  showed  a  method  for 
computer  processing  a  three-dimensional  structure  for  dis- 
play.  A  thesis  by  J.  H.  Kahrs  [Ref.  4]  applied  Desens' 
work  to  displaying  the  picture  seen  by  the  pilot  [Fig.  3] . 

Figure  5  shows  an  object  axis  (aircraft  axis) ,  an 
inertial  axis  (pilot's  eye  glide  slope  intersection  with 
carrier's  centerline) ,  and  a  viewing  plane  axis  (the  LSO's 
eye) .   The  following  set  of  matrices  were  used  to  transfer 
a  set  of  points  described  in  the  object  axis  to  a  set  of 
points  in  the  viewing  plane  axis.   The  matrices  will  also 
give  the  view  a  three-dimensional  perspective  and  control 
its  magnification  (scale) . 

•       [V  V  V  V  = 


[XOB'YOB'ZOB'1] 


Object 
Axis 
Rota- 
tion 


Object 

View- 

View- 

Axis 

ing 

ing 

Trans- 

Plane 

Plane 

lation 

Axis 

Axis 

Trans . 

Rota- 

—          m 

.tion  J 

Perspec- 
tive 
Offset 
Scale 


To  obtain  the  display  coordinates,  division  of  the 
homogeneous  transformed  coordinates  must  be  made  by  the 


W„  coordinate .   Thus 
ri 


Y  =    Y    /W 

DISPLAY  W     H 


ZDISPLAY    :      ZH//WH 
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A. 


Where 

1.  Object  Axis  Rotation  is 
A1  =  cos  3  cos  a  B  p    R 
A2  =  cos  3  sin  a 

A,  =  -  sin  3 

B,  =  -  cos  y  sin  a  +  sin  y 

sin  3  cos  a 

B~  =  cos  a  cos  y  +  sin  Y  sin  3  sin  a 
B_  =  sin  y  cos  3 

C,  =  sin  y  sin  a  +  cos  y  sin  3  cos  a 
C2  =  -  sin  y  cos  a  +  cos  y   sin  3  sin  a 
C,  =  cos  3  cos  y 

2.  Object  Axis  Translation 
10     0     0 
0     10     0 
0     0     10 


XOB   Y0B   Z0B   1 


0 
0 
0 

1 


3.   Viewing  Plane  Translation 
1    0    0    0 
0     10     0 
0     0     10 


"xvp  yvp  ~zvp  1 


4.   Viewing  Plane  Rotation 

Same  as  object  axis  rotation  with 

A,  =  cos  a  cos  3 

A?  =  -  sin  a  cos  y  +  sin  Y  cos  a  sin  3 
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A_  =  sin  a  sin  y  +  cos  y   sin  3  cos  a 

B,  =  cos  B  sin  a 

B?  =  cos  y  cos  a  +  sin  y   sin  B  sin  a 

B.,  =  -  sin  y  cos  a  +  cos  y   sin  B  sin  a 

C,  =  -  sin  B 

C?  =  sin  Y  cos  B 

C~  =  cos  y  cos  B 


5.   Perspective-Offset-Scale 

1     -Y  /F     -Z  /F     -S/F 
o'        o        ' 

0      1        0         0 

0      0         1        0 

0      Y         Z         S 
o       o 

Using  offset  (Y  ,  Z  )  it  is  possible  to  enlarge 

any  section  of  the  picture.   This  feature  is  not  used 

and  therefore 

Y  =  0 
o 

Z   =  0 
o 

The  scale  factor  (S)  used  for  magnification  was 
set  at  0.5  (full  scale).   The  focal  length  (F) ,  used  for 
proper  perspective  was  set  for  the  distance  between  the 
viewer  and  the  surface  of  the  CRT  scope.   The  field  of 
view  is  equal  to 

tan"1 (S/F) 

The  translation  and  rotation  matrices  were  applied 
to  the  problem  of  displaying  the  aircraft  on  the  LSO ' s 
scope  as  follows: 


14 


a.  The  object  axis  was  the  axis  system  containing 

the  lines  necessary  to  draw  the  aircraft.   See  Fig.  7  for 

the  XD,  YD,  Z^  axis.   See  Appendix  F  for  the  proper  method 
a       a       a 

of  describing  the  outline  of  the  aircraft  on  data  cards. 

b.  These  points  were  rotated  to  their  new  posi- 
tion using  the  object  axis  rotation  matrix  with: 

a  =  -ij> 

3  =  a   +  6   +  6 
o    o 

Y  =  <J> 
Since  the  system  was  designed  to  analyze  the 
LSO ' s  function  in  the  loop,  it  was  assumed  that  the  pilot 
will  only  make  small  corrections  to  his  initial  configura- 
tion and  therefore  small  angle  approximations  were  used 
to  cut  down  on  loop  time . 

c.  The  object  axis  was  then  translated  along  the 
vector  OA  shown  in  Fig.  7. 

d.  The  eye  of  the  LSO  (viewing  plane)  was  then 
translated  to  the  LSO  platform.   See  Fig.  7. 

e.  The  viewing  plane  was  then  rotated  so  the  LSO 
is  permitted  to  view  the  approach  from  any  angle  (a) 
[Fig.  5].   ALPHA  equal  to  zero  corresponds  to  looking 
directly  aft. 

Therefore,  in  the  viewing  plane  rotation: 

B  =  0 

Y  =  0 

a  =  YVPLSO 
This  angle  of  rotation  (YVPLSO)  is  fixed  for 
each  approach. 
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The  view  of  the  carrier  deck  for  the  pilot's  display 
is  similarly  generated  where: 

a.  The  axis  describing  the  carrier  deck  is  the 
object  axis  (X  ,  Y  ,  Z  )  see  Fig.  7. 

b.  The  object  axis  was  not  rotated.   Therefore, 
a  =  0,  3  =  0,  y  =  0- 

c.  The  object  axis  translation  was  the  vector  AO. 
(See  Fig.  7)  ■ 

d.  The  viewing  plane  was  not  translated.   Therefore, 

XVP  =  YVP  =  ZVP  =  °' 

e.  The  viewing  plane  is  rotated  where: 

a   =   -    ty  B  =  5        Y  -  ~  <J> 

Any  line  that  extends  behind  the  viewer's  eye  is 
chopped  at  its  intersection  with  the  plane  of  the  CRT 
screen.   Lines  that  extend  outside  the  square  on  the  pilot's 
screen  [Fig.  3]  or  extend  off  the  scope  on  the  LSO's  screen 
are  processed  through  a  software  window  chop  [Ref.  3]. 

Once  the  lines  to  be  displayed  were  processed  through 
the  matrices,  special  software  subroutines  packed  the  dis- 
play buffers  of  the  graphics  computer  and  gave  the  signal 
to  update  the  display.   Figure  6  shows  the  various  steps 
the  digital  computer  performs  in  converting  the  outputs  of 
the  analog  into  pictures  displayed  to  the  pilot  or  LSO . 

C.   ANGLE  OF  ATTACK  INDICATOR 

The  angle  of  attack  indicator  and  the"meatball" [Fig .  3] 
were  displayed  as  in  J.  H.  Kahr ' s  thesis  [Ref.  5] .   The  five 
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states  of  the  angle  of  attack  indexer  are  shown  in  Fig.  6. 

V 


/\     /\ 

1.  _.       2.         3.         4.         5. 


State 
1 
2 
3 
4 
5 


Meaning 
Slow 

Little  slow 
Optimum 
Little  fast 
Fast 


Range  of  AOA 
Greater  than  2.0° 
1.0°  to  2.0° 
1.0°  to  -1.0° 
-1.0°  to  -2.0° 
Less  than  2.0° 


Figure  8.   Angle  of  Attack  Indexer  States 

D .   MEATBALL 

The  meatball  was  visible  within  1.5  degrees  of  the 
optimum  glide  slope.   The  meatball  presentation  [Fig.  3] 
was  located  in  the  middle  of  the  screen  for  an  easier  pilot 
scan.   Only  the  center  square  moved  on  the  screen.   The 
rectangular  shapes  on  either  side  of  the  meatball  repre- 
sented datum  lights  and  did  not  move. 
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III.   DISCUSSION 

The  primary  obstacle  in  modeling  the  system  was  the 
time  required  to  regenerate  an  updated  display.   Several 
time-saving  features  were  used  to  keep  the  regeneration 
time  less  than  .065  seconds. 

1.  A  fast  subroutine  to  pack  the  display  buffers  was 
used  (RTOI) . 

2 .  Small  angle  approximations  were  used  to  draw  the 
aircraft. 

3.  The  number  of  lines  used  to  draw  the  aircraft  was 
kept  to  a  minimum. 

4.  The  flight  dynamics  of  the  aircraft  was  solved  on 
an  analog  computer . 

Hidden-line  removal  [Ref.  3]  would  have  been  an  aid 
for  drawing  a  realistic  picture  of  the  aircraft;  however, 
its  incorporation  would  have  increased  the  loop  time  by 
several  milliseconds. 

Two  expensive  hardware  update  packages  for  the  graphics 
computers  would  have  diminished  the  loop  time  considerably 
so  that  other  options  could  be  included  in  the  program. 
The  hardware  package  to  window  the  drawing  and  the  hard- 
ware to  three-dimensionally  rotate  an  axis  system  would 
have  been  helpful  since  the  software  equivalence  of  these 
items  comprised  a  major  portion  of  the  loop  time. 
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Several  potentially  enhancing  additions  were  not  in- 
cluded in  this  feasibility  study.   An  audible  noise 
generator  linked  to  the  pilot's  throttle  would  aid  the 
LSO  in  determining  the  pilot's  response  to  his  instruc- 
tions.  A  visual  representation  of  the  aircraft's  exhaust 
would  be  helpful  for  the  same  reason.   Isolating  the 
pilot  and  LSO  stations  would  decrease  distractions. 

Once  an  approach  is  completed,  the  LSO  critiques  the 
pass  and  his  judgment  can  be  checked  against  graphs  of 
the  actual  glide  slope  or  lineup  which  are  displaced  on 
the  scope.   This  feature  enables  one  to  determine  the 
ability  of  the  LSO  to  judge  correctly  the  actual  approach, 
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IV.   RESULTS  AND  CONCLUSIONS 

The  LSO  display  was  created  and  tested  to  see  if  devia- 
tions from  the  optimum  glide-path  could  be  detected.   Each 
approach  was  run  on  a  constant  glide-path  which  could  be 
above,  below,  right,  or  left  of  optimum.   The  aircraft's 
angle  of  attack  could  also  be  changed  for  each  pass. 

This  preliminary  model  was  tested  using  five  subjects 
who  were  not  LSO's.   Each  subject  was  shown  several  optimum 
approaches  and  several  approaches  with  non-optimum  condi- 
tions.  The  student  was  then  given  approaches  in  which  he 
had  to  recognize  the  non-optimal  condition. 

These  tests  showed  that  the  subjects  could  judge  lineup 
within  20  feet  of  center  line.   The  difference  between  a 
ramp  strike  and  a  bolter  was  easily  recognized.   Even  angle 
of  attack  changes  of  five  degrees  or  more  were  noticed. 
Although  the  subject  was  unable  to  discern  deviations  from 
optimal  for  every  approach,  his  relative  difficulty  with 
three  parameters  (lineup,  air  speed,  glide-slope)  was  the 
same  as  the  results  of  a  survey  questionnaire  given  to 
LSO's  which  ranked  these  parameters  from  least  to  most  dif- 
ficult to  recognize  [Ref.  5],   Deviations  in  altitude  were 
harder  to  identify  than  for  an  actual  approach. 

The  large  threshhold  of  distinction  for  parameter 
changes  and  the  lack  of  ability  to  determine  altitude  was 
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attributed  to  the  lack  of  training  of  the  student  and  the 
lack  of  a  dynamic  aircraft. 

The  analog  solution  of  the  flight  dynamics  of  an  A-7E 
was  used  to  give  the  aircraft  the  needed  motion.   The  abili- 
ty of  the  analog  computer  to  accurately  describe  a  specific 
aircraft  was  not  germane;  however,  the  model  was  checked 
for  proper  correlation  with  the  A-7E.   Correlation  of  long 
and  short  term  modes  and  the  ratio  of  an  aileron  step  input 
to  the  change  in  velocity  in  the  Z '  direction  are  shown 
below. 

ACTUAL       ANALOG  MODEL 
Long  term      65.8  sec      65.1  sec 
Short  term     4.85  sec      4.85  sec 
6e/DW  3.17  3.0 

The  interaction  among  the  pilot,  the  chair,  and  the  pre- 
sentation on  the  pilot's  scope  was  tested  by  the  J.  H. 
Kahr's  [Ref.  4].   Kahr ' s  concluded  that  "sufficient  visual 
cues  are  provided  to  enable  consistent  landings  by  experi- 
enced pilots . " 

The  simulation  presented  in  this  thesis  demonstrated 
the  feasibility  of  using  a  hybrid  computer  linked  to  graphic 
output  to  present  an  environment  where  the  interaction  of 
the  LSO  in  the  carrier  approach  landing  system  could  be 
studied.   This  program  could  be  used  as  a  basic  tool  to 
study  the  LSO/pilot  interaction;  however,  improved  hardware 
capabilities  would  enable  more  extensive  studies. 
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Figure  2.   Allocation  of  Computer  Functions 
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Figure  3.   Pilot  Display. 
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Figure    4.      LSO   Display, 
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APPENDIX  A 
SUMMARY  OF  EQUATIONS  OF  MOTION 

Etkin's  [Ref.  2]  longitudinal  and  lateral  equations 
are  reprinted  below: 

Longitudinal 

(2yD  -  2CT   tan  9^  -  C   )u  -  C   a  +  CT   6  =  0       (1) 
Lo      °    Xu      xa      Lo 

(2C    +  C   )u  +  (2yD  -  C    -  C   )a 
o     u  a     a 

-[(2y    +   Cz    )D   -   CL      tan    9q]9    -   C%   r\   =   0  (2) 

q  o  n 

-C        u   -     (C        D   +   C      )o   +    (i_    D2    -    C         D)9 

m  m*  m  B  m 

u  a  a  q 

-    (C        D  +  C      )    n   =   0  (3) 

n  r\ 


where    V     ps£ 

*         2. 

1  -k 

%        2 

Lateral 

(2yD   -   C      )B    -   C        p   +    (2y    -   C      )r   -   C        * 
y3  yp  yr  o 


+  C        ?   =   0  (4) 


-    %    6    +     (iA   D    "    C£    }P    "     {iE   D   +   C*r)r 


-    (C,       D    +    Cp     K    -    C.    5    -    0  (5) 

J6?  *?  a 


-C       B  -(iE  D  +  cn   )P  +    (ic  d  -  cn   )r 
@  P  r 


-C         g    -    (C         D    +   C)^    =    0  (6) 

n5  n^  nc 
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The  following  substitutions  and  assumptions  were  used 
to  reduce  these  equations  to  a  more  usable  form. 
In  General 


*   d 

D  =  t   at 

u 

u                               u 
o                              o 

U   =   u   +   u 
_o 

P ,      =   i  p   u 
d^        2    r      o 
o 

Longitudinal 

I      (y) 

i    -    yy  Y 

b             c   -3 
p    S   c 

2m                       * 

y   =  — -                   t     = 

pSc 

c 

2u 
o 

Equation   1 

(3T/8uK 

c       -                 ° 

xu           Pd      S 
o 

-    2  (C_      +   CT      tan   6    ) 
o            o 

-   M 

Cx      -   CL      "    C0 
a            o            a 

9M 


where 


OT/3u)  3Cn 

P,    S       =  °      and     M  5iT  =  ° 

d, 


0 

Equation  2 

3C 
czu  -  -  M  8FT  =  ° 

CZ   "  "(CL   +  CD  > 
a       a     o 


cz-  =  °      cz  -  ° 

a  q 
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Equation  3 

\ 

=  0 

=  c 

6e 

Lateral 

P  = 

2u 
o 

qc 

rb 

r  "  2u 
o 

i  = 
a 

I  (y) 

XX  ' 

p  Sb3 

zz 

lc    p  SB3 

Lc  = 

I  (y) 

xz  ' 
p  Sb3 

2m 

y  "  PSb 

2u 
o 

Equation  4 

C 

X,  -   c 

y6r 

6r 

Equation  5 

C^c 

=  C^ 

CJU  =   ° 

I 

6a 

K 

5 

fir 

Equation  6 

- 

% 

=  C 

n. 
6a 

.    «*•  .  0 

n     n .. 
p     6r 

As  a  result  of  the  above  substitutions  and  assumptions, 
the  following  equations  resulted  and  were  solved  by  an 
analog  computer. 
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The  remaining  equations  are  reprinted  from  p.  122  Ref.  2. 
6  =  q  (13) 

<f>  =  p  +  r  tan  6  (14) 

$  =  r  sec  6Q  (15) 

o 

x  =  u   cos  6   +  u  cos  6   -  u  6  sin  9   +  6   w        (16) 
o      o  o    o       o    o 

o 

y'=  u  i>   cos  0   +  v  (17) 

2  o  r      o 

o 

z'=  -  u   sin  6   -  sin  0  u  -  u  0  cos  0   +  w  cos  0    (18) 
o      o        o     o       o  o 

Equations  7-18  were  scaled  and  solved  on  the  Cl-5000 

analog  computer  as  shown  in  Appendix  B. 
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APPENDIX  B.   THE  ANALOG  PROGRAM 
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A052 


A003 
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P055 
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A007 
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T505 


P413 
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T500 
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T502 


Z' 


36 


T121 


P057 


Thrus 


P405 


P027 
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A206 


A064 


P040 
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fir 


T104 
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P,   S  b  C„   SCV 
a       x,  0 

P001  =    °        '; 


P002  = 


I 

XX 

u 
o 

SCP 

pd 

o 

P 

Sb2 

2u 
o 

I 

XX 

(10) 

• 

pd 

o 

*b2 

C£   SCR 

r 

I 

XX 

2u 
o 

SCP 

I 

XX 

SCP 

I 
xz 

SCR 

I 

XX 

SCP 

P003 


P,   S  B  C0    SCDA 
a        x-j, 

XX 


P,   S  b  C0    SCDR 

d  £  r. 

POOS  =    °       6r 


P006  = 


P ,   S  b  C    SCV 
d^       nR 

P007  =  — 2 E 

I      u      SCR 
zz    o 

P,   S  b2  C    SCP 
d        n 

P010  =  — E 

ruxu    I    2u        SCR 
zz    o 

P,   S  b2  C 

d        n 

POii  =  _° L 

I    2u 
zz    o 


P ,   S  b  C     SCDA 
d       n- 

P012  =  —2 iS 

*u      I  SCR 

zz 


P ,   S  b  C     SCDR 
d       n  r. 

P013  =        °       6r 


Izz  SCR(IO) 
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I 
xz 

SCP 

I 
zz 

SCR 

pd 

o 

S  C 

m 

u 
o 

Pd 
o 

S  b  C 

YP 

SCP 

2 

m  u 
o 

scv 

P014  = 


P015  = 


P016  = 


P ,   S  b  C    SCR 
d       y 

P017  =  —2 1* 

2    m  u      SCV 
o 


P,   S  C     SCDA 

p°2°  ■  TT ^scv- 


P ,   S  C     SCDR 
d      Vj 

P021  =    °      6r 

*^X    m      scv  (10) 


P-   S  CT   SCPH 
d      L 

p022  =  ~ST °scv(io) 


-  u   SCR 

P023  =  scv°(io) 


Pd   S  (-2)  CD 

P024  =  —2 ° 

m  u 
o 


Pd   S  (CL   "  CD  }  SCW 

P025  =  — - ^ 2 

SCU 


P026  = 


P027 

m  SCU 


m 

u 
o 

"Pd 
o 

S  CL 
o 

SCTH 

m 

SCU 

SCDT 
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-2  P,   S  CT   SCU 
a     L 

P030  =  2 °^ 

m  u         SCW 
o 


"  <CL   +  CD  >  Pd   S 

P031  =  2 2 2_ 

m  u 
o 


-  CT-   P,   S  (tan  6  )  SCTH 

J_i     Q  O 

P032  =  2 2 

m  SCW 


uQ  SCQ 
P033  =  SCW (10) 


P ,   S  c  C    SCU 
d_      m 

P034  =  —^ H 

I    u        SCQ 
yy   o 


P ,   S  c  C   SCW 
d       m 

P035  =  —2 2 

1  u       SCQ 
yy   o 

P,   S  c2  C    SCW 
d        m* 

P036  =  —2 ^—2 

2  Jyy  uo     SCQ 

P,   S  c2  C 

d        m 

P037  =  x-s a 

2  1    u 
yy   o 


P  ,   ScC     SCDE 
P040  =   °       6e 


I  „  SCQ  (10) 

yy 


tan  8   SCR 


P041  = 


SCPH 


sec  6   SCR 

P042  = 

SCSI 


u   cos  0 
P413  =  -2 2 

sex 
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cos  e 

o 

SCU 

sex 

-  u  sin  0 
o      o 

SCTH 

P044  = 


P045  =       SCX 


e  sew 
p046  =  "lex— 


u   cos  6   SCSI 
P047  =  -2- 


SCY 

-  u   sin  6 

P417  =  2 ° 

SCZ 

P050  =  .50 


P051 


P052  = 


-  sin  8 
o 

SCU 

SCZ 

-  u  cos 
o 

6 
o 

SCTH 

SCZ 

cos  6   SCW 
o 

P053 

SCY 

P054 

= 

SCQ/SCTH 

P055 

= 

SCP 

(10)SCPH  . 

P056 

= 

SCU/SCY 

P057 

= 

.50 

P405 

= 

.1716  (BIAS) 

P406 

= 

-.0576 (BIAS) 

P407 

= 

-.0963 (BIAS) 

P416 

= 

-.0594 (BIAS) 
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APPENDIX  C.   DIGITAL  PROGRAM  FLOW  CHART 
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APPENDIX  D 
PROGRAM  VARIABLES 


AIXX 

AIXZ 

AIYY 

AIZZ 

ALPHA  <j> 

AMBX 

ANGL 

AOA 

B 

CB 

CDA 

CDO 

CHECK 

CLA 

CLO 

CNST 

1 

CNST 

2 

CNST 

3 

CNST 

4 

CNST 

5 

-  Moment  of  inertia  about  the  XD  axis 

-  Product  of  inertia 

-  Moment  of  inertia  about  the  YD  axis 

-  Moment  of  inertia  about  the  Z  axis 

-  Initial  angle  of  attack  (a  ) 

-  Array  of  vertical  position  of  meatball 
lines 

-  Actual  glide  slope  angle  of  the  aircraft 

-  Angle  of  attack 

-  Wing  span 

-  Mean  aerodynamic  chord  (c) 

-  C^ 


a 


-  c. 


X'  distance  from  LSO  to  the  starting 
position  of  the  aircraft 


-  C, 


-  C, 


a 


CTHO 


"  a0  +  eo 

-  CNST  1+9 

-  Maximum  position  of  the  meatball  before  it 
goes  off  the  lens 

-  0.5/75.0 

-  Distance  from  starting  point  of  the  A/C 
to  the  stern  of  the  carrier 


cos  (8  ) 
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-  cos  (Y) 

-  Difference  between  actual  glide  slope 
angle  and  ideal  glide  slope  angle 

-  Focal  length  for  use  in  prospective  matrix 

-  Flag  which  is  set  when  aircraft  is  past 
.  the  LSO 

-  Horizontal  intersection  of  a  line,  used 
in  the  software  window 

-  Integer  variable,  usually  a  counter 

-  Number  read  from  LSO ' s  input  to  determine 
whether  the  restart,  present  a  graph,  or 
terminate  option  is  wanted 

-  First  point  in  XKEEP  array  with  a  value 
less  than  2000  ft 

-  Number  of  data  cards  needed  to  draw  the 
aircraft 

-  A  counter  used  to  check  the  number  of 
times  a  point  has  gone  through  the  software 
window 

-  A  packed,  bogus  array  used  to  maintain  con- 
stant loop  time  after  the  aircraft  passed 
the  LSO 

-  Packed  array  for  a  fast  AOA  display 

-  Graphics  directory  for  AGT  #n.   Can  use 
eight  graphics  blocks  maximum 

-  Packed  array  for  a  slow  AOA  display 

-  Move  draw  array  for  the  meatball 

-  Packed  array  for  an  on  speed  AOA  display 

-  Packed  array  with  the  fixed  display  shown 
to  the  LSO 

-  Array  of  fixed  coordinates  for  graphical 
glide  slope  or  lineup  display 
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IPL0T2 

IPL0T3 
ISQ 

ITDIRn 

ITEST 

ITEXnn 

IX/IY 


L 
LL 
M 
MDC 

NAC 

NB 

NL1 

NL2 

NL3 

NL4 

NLAND 

NLC 


-  Packed  array  to  store  actual  glide  slope 
track  for  graphing 

-  Packed  array  to  store  actual  lineup  track 

-  Packed  array  for  fixed  graphics  shown  to 
the  pilot 

-  Text  directory  for  AGT  #n.   Maximum  of  19 
text  blocks  allowed 

-  Decision  variable  to  read  in  new  data 
information  or  terminate  the  problem 

-  Line  of  text  sent  from  the  digital  computer 
to  the  graphics  computer 

-  Operating  values  of  start/end  points  used 
in  software  window 

-  Number  of  points  to  be  plotted  on  the 
graphs 

-  Integer  variable  used  as  a  counter 

-  Counts  once  for  every  10  regenerating  loops 

-  Aircraft  mass 

-  Move  draw  array  for  the  runway  and  meat- 
ball presentation  to  the  pilot 

-  Number  of  data  cards  necessary  to  describe 
the  aircraft 

-  2*  NL4  +  1 

-  NLSHIP  +  1 

-  NLSHIP  +  2 

-  NLSHIP  +  3 

-  NLSHIP  +  4 

-  Number  of  lines  needed  to  draw  landing  area 
as  seen  by  the  pilot.   One  line  per  data 
card 

-  NLAND  +  1 
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NLC1 

-  2  *  NLC 

NLC2 

-  2  *  NLC  +  1 

NLC5 

-  NLC1  +  5 

NLC6 

-  NLC1  +  6 

NLSHIP 

-  Number  of  d 

ONEDEG 

ONEHAF 
P 

PHI 

PSI 

Q 
S 

SCDA 
SCDE 
SCDR 
SCDT 
SCTH/SCPH/SCSI- 

SCU/SCV/SCW 

SCX/SCY/SCZ 

SLOPE 


Nuniber  of  data  cards  needed  to  represent 
the  lines  which  outline  the  carrier  deck 
as  seen  by  the  LSO 

-  Limit  for  AOA  display.   Equals  one  degree 
in  radians 

-  Limit  for  meatball  display 

-  Working  matrix  for  aircraft  stability 
derivatives 

-  Scaled  <f>  (change  in  roll  angle  from 
initial  conditions) 

-  Scaled  ij>  (change  in  yaw  angle  from 
initial  conditions) 

-  Dynamic  pressure 

-  Wind  area  of  the  aircraft.   Also  used  as 
the  scale  factor  in  the  prospective  matrix 

-  Scale  factor  for  maximum  aileron  deflection 

-  Scale  factor  for  maximum  elevator  deflection 

-  Scale  factor  for  maximum  rudder  deflection 

-  Scale  factor  for  maximum  thrust 

Scale  factors  for  maximum  values  of 
6/<}>/and  \\> 

Scale  factors  for  maximum  velocity  change 
in  the  u/v/w  directions 

Scale  factors  for  maximum  travel  in  the 
X'/Y'/Z1  directions 

Slope  of  display  line,  used  in  software 
window 
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STERN 

SYAW 
STHO 
TEMP1+2 

THETA 

TH<j> 

TWODEG 

U<(> 
VINSCT 

W 

WIREn 
X/Y/Z 
X0/Y0/Z0 


Distance  from  the   FRESNEL  LENS  to  the 
Round-down 

sin  (Y) 

sin  (TH<J>) 

Temporary  storage  for  a  point  on  the 
horizon 

6  (change  in  angle  of  climb  from  the 
initial  condition) 

6, (initial  angle  of  climb) 

Limit  for  AOA  display.   Equals  two 
degrees  in  radians 

u   (initial  aircraft  velocity) 

Vertical  intersection  of  a  line  used 
in  the  software  window 

Scaled  change  of  velocity  in  the  Z ' 
direction 

Distance  from  the  FRESNEL  LENS  to  the 
nth  wire  x,  n  =  1,2,3,4 

Scaled  inertial  X'/Y'/Z1  positions  of 
the  aircraft  from  its  starting  position 

Position  on  the  carrier  deck  where  the 

LSO ' s  eye  is  located  (measured- in  the 

X  ,  Y  ,  Z   axis  system  [Fig.  5]) 
c   c   c       J         - 


XAC/YAC/ZAC    -  Coordinate  system  in  which  the  aircraft 

is  drawn  (X^.Y^/Z   in  Fig.  5) 

D    D    D 

XACn/YACn/ZACn-  Working  vectors  for  aircraft  data 

XEn/YEn/ZEn   -  Working  vectors  to  handle  carrier  deck 

points 


XELND/YELND/ 
ZELND 


End  point  for  lines  that  represent  the 
landing  area  as  seen  by  the  pilot 


XESHIP/YESHIP/ 
ZESHIP 


Working  vectors  used  to  process  carrier 
deck  graphics  as  seen  by  the  LSO 


51 


XESHP/YESHP/ 
ZESHP 


End  point  of  lines  describing  the  carrier 
deck 


XINIT/ZINIT   -  Starting  position  for  the  aircraft 

(measured  in  X  ,Y  ,  Z   axis  system  with 
Y   =  0.0)      c   c   c 
c 


XKEEP/YKEEP/ 
ZKEEP 


XSn/YSn/ZSn 


XSLND/YSLND/ 
ZSLND 


XSSHIP/YSSHIP/ 
ZSSHIP 


XSSHP/YSSHP/ 
ZSSHP 


XT/YT/ZT 
XTl/YTl/ZTl 


YMB 

YTEMP 

YUPLSO 
Z 

ZINIT1 


Array  holding  inertial  X/Y/Z  position  of 
the  aircraft  for  use  in  the  glide  slope  and 
lineup  graphs 

Working  vectors  used  to  handle  carrier 
deck  points 


Starting  point  for  lines  used  to  draw  the 
landing  area  as  seen  by  the  pilot 


Working  arrays  used  to  process  carrier  deck 
graphics  as  seen  by  the  LSO 


-  Starting  point  of  a  line  used  to  describe 
the  carrier  deck  as  seen  by  the  LSO 

-  Inertial  distance  from  the  LSO  to  the 
aircraft  (X'^Y'^Z1) 

-  Inertial  distance  from  the  pilot  to  the 
origin  of  the  X  ,Y  ,Z   axis  [Fig.  5] 

\*r  W        \-f 

-  Value  used  in  software  window,  also  used 
as  the  scaled  Y1  position  of  the  aircraft 

-  Array  containing  the  horizontal  position  of 
the  displayed  meatball 

-  Temporary  storage  for  Y-coordinate  in 
software  window 

-  Yaw  angle  of  the  viewing  plane  of  the  LSO 

-  Value  used  in  software  window,  also  the 
scaled  inertial  Z'  position  of  the  aircraft 

-  Distance  from  aircraft  starting  point  to 
the  water 


ZSAVE 


Save  the  value  of  Z  because  of  its  dual 
use  in  the  window 
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APPENDIX  E 
OPERATING  INSTRUCTIONS 

The  program  exists  in  two  forms.   A  card  deck  which  must 
be  input  through -the  card  reader  is  available  in  the  computer 
room  in  a  box  marked  "Pilot/LSO/aircraft  interaction."   The 
program  is  also  located  on  a  tape  marked  in  the  same  manner 
as  the  box.   This  tape  is  a  core  dump  and  doesn't  give  a 
listing  when  used.   Use  of  the  tape  requires  about  15-20 
minutes  less  to  set  up. 

Operating  instructions  for  the  various  computers  is  lo- 
cated in  the  computer  room  on  the  fifth  floor  of  Spanagel 
Hall. 

1.  Load  the  tape  marked  "Pilot/LSO/aircraft  interaction" 
or  place  the  large  deck  of  cards  in  the  card  reader. 

2.  On  the  Tape  Control  Panel:   Put  tape  at  the  load 
point,  select  automatic,  set  unit  select  to  #1. 

3.  If  the  large  card  deck  is  used,  mount  a  scratch  tape 
and  omit  step  4. 

4.  Put  the  card  deck  marked  "For  use  with  tape"  in 
the  card  reader.   Ready  the  card  reader. 

5.  Push  IDLE,  RESET,  RUN,  CARDS  on  the  XDS  9  300  con- 
trol panel. 

6.  Mount  analog  and  logic  boards  numbered  21  onto  the 
CI  5000. 
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7.  Press  "Digital  Computer"  on  the  analog  control 
panel.   The  only  lights  that  should  be  on  are  LOCAL,  POTSET, 
IDC-X1,  REAL  TIME,  and  CL0CK1-H. 

8.  Seven  pots  must  be  set  manually.   A  list  of  these 
pots  and  their  values  should  be  listed  on  the  XDS  9  30  0 
teletype.   Set  the  seven  pots.   Return  the  analog  to  "Digi- 
tal Control . " 

9.  Set  the  pilot's  chair  in  front  of  the  AGT  #2. 

10.  Ready  both  AGT  units  and  set  in  GATDl  by  typing: 

RESET  ("GATDl",  104) ! 
GATED! 
If  "FILE  NOT  FOUND"  appears,  type  04  vice  104. 

11.  The  program  asks  the  operator  to  type  a  one  (1) 
on  the  XDS  9300  teletype  "to  continue."   The  program  will 
have  to  be  reloaded  if  both  AGT  units  are  not  ready  when  a 
1  is  typed.   Type  a  1  and  a  carriage  return  on  the  XDS  9  300 
teletype . 

12.  Text  should  appear  on  both  AGT  scopes.   At  this 
point  the  AGT-1  teletype  has  control  of  the  problem.   Follow 
the  instructions  on  the  scope. 

13.  After  an  approach  has  been  completed,  the  LSO  has 
three  options : 

a.  Begin  a  new  approach. 

b.  Have  displayed  a  graphical  output  of  glide  scope 
and  lineup  results  from  the  approach. 

c.  Terminate  (i.e.  give  control  back  to  the  proctor) 
If  terminate  is  selected,  the  following  will  appear  on  the 
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XD5-9  300  teletype:   "Type  1  to  Read  new  data,  Type  0  to 
Stop".   The  proctor,  at  this  point,  may  change  any  data 
card  and  place  the  data  deck  (without  the  control  cards) 
in  the  card  reader.   Then  type  a  1  on  the  XDS-9  300  tele- 
type and  go  to  step  11. 

14.   The  -program  is  ended  by  typing  a  zero  instead  of 
a  one  on  the  XDS  9  30  0  teletype. 
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APPENDIX  F 
PREPARATION  OF  DATA  DECK 

Card  1-7 

The  first  seven  cards  are  used  to  read  in  non-dimen- 
sional aircraft  stability  derivatives  and  other  aircraft 
parameters.   Ten  spaces  are  provided  for  each  number. 
Slots  1-10  are  used  for  the  first  number  on  a  card,  11-20 
for  the  second  number,  etc.   The  first  seven  cards  and 
their  contents  are  listed  below.   Note:   All  non-dimensional 
derivatives  are  per  radian  except  as  noted. 

1.  5 (ft),  6 (ft),  c(ft),  I   (slug-ft2) ,  I   (slug-ft2) , 

xx  yy  , 

I   (slug-ft2),  I   (slug-ft2) 
zz     3        xz 

2.  Cg    ,  C0  ,  C«  ,  C.    (per  deg) ,  Cp    (per  deg) 

*g    p    r   *6a  *6r 

3.  C   ,  C   ,  C   ,  C     (per  deg),  C     (per  deg) 

ll/i       11         11        11  c  11  I" 

3    p    r    6a  or 

4.  C   ,  C   ,  C   ,  C     (per  deg),  C     (per  deg) 

yB   yp   yr   y6a  y6r 

5.  C   ,  C   ,  C   ,  C   ,  C     (per  deg) 

m    m  '   m*    m  '   m~    r     s 
u    a     a    q     oe 

6'   CL  '  CD  '  CD  '  CL  '    Uo  (ft/Sec) 
a    o    a    o 

7.   6  ,  A/C  mass  (slugs) 

Note:   Use  a  decimal  point  with  all  numbers. 

Card  8&9 

The  next  two  cards  contain  scaling  values.   Again,  ten 

spaces  are  provided  for  each  number.   Use  scale  values  for 
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the  following  variables: 

8.  fir,  Se,  6a,  5    ,    u,  v,  w,  9 

9.  <j>,  ij;,  X'  ,  Y'  ,  Z'  ,  p,  q,  r 

10.  ALPHA  zero  in  the  first  10  spaces. 

11.  The  next  set  of  cards  describe  the  lines  which  repre- 
sent the  carrier  deck  seen  by  the  LSO.  The  first  card 
contains  the  number  of  lines  needed  to  draw  the  ship 
(NLSHIP-I2) .   NLSHIP  is  the  number  of  cards  that  follow 
since  each  line  used  to  draw  the  ship  uses  one  card. 

The  next  NLSHIP  cards  are  prepared  using  the  X  ,  Y  ,  Z 

c     c  a      c   c   c 

axis  system  [Fig.  5]  with  its  origin  fixed  to  the  carrier 
deck  as  the  point  where  the  center  line  of  the  angled  deck 
meets  a  line  perpendicular  to  the  center  line  and  going 
through  the  FRESNEL  LENS.   The  positive  X   axis  points 
toward  the  bow  of  the  ship,  the  positive  Y   axis  points 
toward  the  port  side  of  the  ship. 

Expressed  in  the  X  ,  Y  ,  Z   axis  system,  the  lines  de- 
scribing the  carrier  deck  are  put  onto  cards  as  follows: 

Columns  1-10,  11-20,  21-30  the  coordinates  of  the 
beginning  of  a  line. 

Columns  31-40,  41-50,  51-60  the  coordinates  of  the  end 
of  the  same  line. 

The  next  set  of  cards  outline  the  aircraft.   The  first 
card  contains  the  number  of  cards  needed  to  describe  the 
aircraft  (NAC) .   Format  12 

The  axis  system  used  is  (Xn,  Y_, ,  Z_)  [Fig.  5]  with  the 

ODD 

origin  at  the  center  of  gravity  of  the  aircraft.   The 
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positive  Xc  axis  points  toward  the  nose  of  the  A/C,  the 

ri 

positive  YD  axis  is  out  the  left  wing  and  the  positive  Z_ 
axis  is  up..  In  this  axis  system,  put  the  coordinates  of 
all  intersecting  lines  in  column  1-10  (X...)  ,  11-20  (Y_.)  , 
and  21-30  (Z  ) .   In  column  31  place  a  zero  to  move  and 
a  1  to  draw 'a  vector. 

Maximum  number  of  lines  is  60. 

The  next  card  contains  the  number  of  cards  following 
which  describe  the  landing  area  as  seen  by  the  pilot. 
Formal  12.   The  program  had  five  lines  drawn  and  represent 
a  rectangular  runway  10  0  0  ft.  long  and  10  0  ft.  wide.   It 
is  recommended  that  these  cards  remain  the  same  since  more 
lines  would  increase  the  loop  time.   If  these  cards  are 
replaced,  the  format  and  axis  system  is  the  same  as  the 
one  used  in  describing  the  carrier  deck  as  seen  by  the  LSO. 

The  next  card  contains  in:  columns  1-10,  11-20,  21-30, 
the  position  of  the  LSO's  eye  expressed  in  the  X  ,  Y  ,  Z 
axis  system;  in  column  31-40  the  scale  factor  (S)  equal  to 
1.0;  in  column  41-50  the  focal  length  (F)  equal  to  20  feet. 

The  next  card  contains  the  number  of  degrees  the  LSO 

moves  his  head  and  this  his  field  of  view.   A  positive  angle 

is  equivalent  to  moving  his  head  to  the  left.   Use  column 

1-10. 

The  next  card:   Column  1-10,  11-20  the  X   and  Y   start- 

'  c      c 

ing  position  of  the  aircraft.  (XINIT  and  YINIT) 
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The  last  card  contains  in  column  1-10,  11-20,  21-30,  and 
31-40  the  distance  between  the  FRESNEL  LENS  and  the  #1,  #2, 
#3,  #4  wire  respectively.   Columns  41-50  contains  the  dis- 
tance from  the  FRESNEL  LENS  to  the  stern  of  the  ship. 
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Wire  1    Wire  2 


Wire  3 


Wire  4 


Stern 


KINIT 


YINIT 


YVPLSO 


XO 


YO 


ZO 


NLAND+NUMBER  OF  CARDS  TO  FOLLOW  WHICH  DESCRIBE 
THE  PILOT  VIEWED  LANDING  AREA 


NAC+NUMBER  OF  CARDS  TO  FOLLOW  WHICH  DESCRIBE 
THE  AIRCRAFT 


11  NLSHIP+NUMBER  OF  CARDS  TO  FOLLOW  WHICH  DESCRIBE 
THE  SHIP  (12) 

10  a 
o 

8&9  SCALING  PARAMETERS  FOR  ANALOG  COMPUTERS 


7.   6  ,  A/C  MASS  (SLUGS). 


1-6  A/C  CONSTANTS 


Data  Deck  Sample 
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APPENDIX  G.   COMPUTER  PROGRAM 
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